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Abstract 

In this review we describe statistical mechanics of quantum systems in the pres- 
ence of a Killing horizon and compare statistical-mechanical and one-loop contribu- 
tions to black hole entropy. Studying these questions was motivated by attempts to 
explain the entropy of black holes as a statistical-mechanical entropy of quantum 
fields propagating near the black hole horizon. We provide an introduction to this 
field of research and review its results. In particular, we discuss the relation between 
the statistical-mechanical entropy of quantum fields and the Bekenstein-Hawking 
entropy in the standard scheme with renormalization of gravitational coupling con- 
stants and in the theories of induced gravity. 
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1 Introduction 



Thermodynamics and statistical mechanics of black holes is one of the most exciting and 
rapidly developing areas of black hole physics. Black holes are known to possess the 
properties similar to the properties of thermodynamical systems According to this 
analogy, a black hole has the entropy S BH = jqA where A is the surface area of its 
horizon and G is the Newton constant^ The quantity S BH was introduced in Refs. 0-0 
and is known as the Bekenstein- Hawking entropy. The temperature of a black hole is 
Th = 7^ where k is the surface gravity of the horizon. 

In the Einstein theory of general relativity the Bekenstein-Hawking entropy is a pure 
geometrical quantity. The laws of black hole thermodynamics are derived by making use of 
the classical Einstein equations and rules of the differential geometry only. If we compare 
black holes with usual thermodynamical systems an important difference can be easily 
observed: Black holes are nothing but an empty space with a strong gravitational field 
while a usual body consists of material constituents (atoms, molecules, etc.). Namely 
this microscopical structure enables one to explain thermodynamical properties of the 
body in terms of statistical mechanics of its constituents. Does a black hole have internal 
degrees of freedom which are responsible for the Bekenstein-Hawking entropy? This is a 
key question of black hole physics. 

The black hole entropy problem is important. The statistical- mechanical derivation 
of the Bekenstein-Hawking entropy is a highly non-trivial test for a fundamental theory 
of quantum gravity. Recent successful calculations of S BH for extremal 0-@ and near 
extremal |10|, [pl]1 black holes in the superstring theory clearly demonstrate this. A review 
of these computations and further references can be found, for instance, in |12|], [ [L3[] . 
Besides the superstring approach there are other approaches which attack the problem of 
black hole entropy from different directions. The explanation of the entropy of 3D black 
holes suggested by Carlip et al. [14|,[15] or consideration in the framework of the loop 
quantum gravity [|l6j are some of them. 

The subject of this review is related to the approach which suggests to explain the 
black hole entropy in terms of its quantum excitations. This idea was formulated first in 
fT7fl,|1l8||,fl9|l and it has stimulated a large number of publications. 

The properties of physical vacuum, especially in the presence of gravity, are nontrivial. 
In the state of vacuum there always exist zero-point fluctuations of physical fields. An 
observer at rest near a horizon would register these zero-point fluctuations in the form of 
a thermal atmosphere of a black hole |19|]-[[25|]. Historically the first suggestions to relate 
the entropy of a black hole to the entropy of its thermal atmosphere were made in works by 
Thorne and Zurek [17| and by t'Hooft [fL8| . t'Hooft |18|] estimated the thermal entropy by 
assuming that the red-shifted temperature of the atmosphere is Tjj and showed that the 

1 We use units in which h = c = fee = 1 is the Bolzmann constant) and sign conventions of book 
[p| and, thus, we use the signature (— , +, +, +) for the Lorentzian metric. 
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entropy is proportional to the horizon area A. To avoid the divergences t'Hooft assumed 
that fields vanish within some distance near the horizon (the corresponding model was 
called the "brick wall" model). If this distance is of the order of a Planck lengths the 
thermal entropy turns out to be comparable to S BH . 

The reason why the static observer near the black hole sees the vacuum as a mixed state 
is explained by the loss of the information about a part of the quantum system located 
inside a the black hole horizon. Bombelli, Koul, Lee and Sorkin [26[ and Srednicki |27| 
have showed that even in a flat space when observations in vacuum are restricted by a 
part of the system, the entropy is not zero and is proportional to the surface area of the 
restricted region Q (for a recent discussion see p8[ ) . Similar result for the entropy was also 
established for non-zero spin fields p9j and for the pure sates different from the vacuum 



30]. The non- vanishing entropy appears because the "observable" and "non-observable" 



vacuum fluctuations are entangled (correlated) on the boundary of Q. By taking into 
account these properties the authors of Refs. fl26|, ]p7H suggested to explain S BH as the 
entropy of entanglement between quantum fluctuations propagating on the different sides 
of the horizon. 

Frolov and Novikov |3T| proposed to relate the black hole entropy to the degrees of 
freedom corresponding to quantum states in the black hole interior. The density matrix of 
these degrees of freedom can be obtained by averaging the quantum state of the complete 
system over the states of fields located outside the black hole. For modes in the close 
vicinity to the horizon this density matrix is thermal. The particles are created in pairs, 
and only one of the components can be created outside the horizon. A pair inside the 
black hole is in a pure state and does not contribute to the entropy. For this reason, the 
statistical mechanical entropy is connected with entanglement, and it can be written in 
the form of summation over the modes in the black hole exterior. In other words, this 



approach incorporates main features of the earlier approaches [17] , [IS ] ,[26] . 

A remarkable property of a black hole is that its entaglement entropy, and entropy 



connected with its thermal atmosphere coincide ||32||-||35||. In what follows we refer to this 
quantity as to statistical-mechanical entropy. 

Small fluctuations of fields (including the gravitational one) propagating in the black 
hole background can be related to small deformations of the black hole geometry. This 
can be explicitly demonstrated in the approach using the no-boundary wavefunction of 
the black hole []36(| . For this reason, counting states of quantum fields is connected with 
the counting the states of quantum excitations of the black hole. 

In the general case, the relation of the statistical-mechanical entropy and the "observ- 
able" thermodynamical entropy of the black hole is highly nontrivial [37]. The quantum 
fields near a black hole have an important property. Namely, the density of levels of quan- 
tum mechanical Hamiltonians of the fields blows out near the horizon. This results in 



the divergence of the statistical- mechanical entropy. Susskind and Uglum [38 and Callan 



and Wilczek [34| suggested that this divergence is related to the ultraviolet one-loop di- 
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vergences of the theory and it can be removed by renormalizing the Newtonian coupling 
constant. This observation was strongly supported by the result of Demers, Lafrance 
and Myers |39| who showed that the standard Pauli-Villars method regularizes the di- 
vergences of the statistical-mechanical entropy on the Reissner-Nordstrom background. 
Moreover in this regularization, the entropy divergence for a minimally coupled scalar field 
is completely eliminated by the standard ultraviolet renormalization. The same result for 



general static black-hole backgrounds was proven by Fursaev and Solodukhin jMj . It was 
done by using the Euclidean formulation of the theory with conical singularities. 

54j is very important, there are two problems. 



Although the observation made in [38 



First, in order to carry out the renormalization one must introduce an infinite (negative) 
bare entropy, which has no statistical-mechanical origin. Second, in the presence of non- 
minimal coupling of a field with the curvature, the divergence of the statistical-mechanical 



entropy and the ultraviolet divergence are different [39] — 144] . 

Jacobson [45[] pointed out that the first problem can be solved if the Einstein gravity 
itself appears as a result of quantum effects. A suitable simple example is Sakharov's 
theory of induced gravity |4(J,[^7|]. The models of induced gravity which enable one to 
check this hypothesis were constructed by Frolov, Fursaev and Zelnikov fl48|. It was 



demonstrated that the presence of the non-minimal coupling is necessary in order to get 
a finite induced Newton constant. Moreover, it makes it possible a consistent derivation 



of the black hole degeneracy by counting the degrees of freedom of constituents [49|,[50 



The idea to relate S BH to quantum excitations of a black hole and the pioneering 



papers Ul7| , flgl , p6| , f27| , fl31[ stimulated the study of statistical mechanics of quantum 
fields in the presence of a Killing horizon. More than hundred papers was written on this 
subject. Our review of these works has two main aims: 

1. To describe the statistical mechanics of quantum fields in the presence of a Killing 
horizon, its methods and obtained results; 

2. To discuss the relation between the statistical mechanical entropy and the observable 
black hole entropy S BH . 

As an application of these results, we present the statistical-mechanical derivation of 
the black hole entropy in the models of induced gravity. 

We begin with two formulations of finite-temperature quantum theory on static space- 
times which have been discussed in the literature. 

The first is the canonical formulation based on the (3 + l)-split of the spacetime. It 
allows one to define the free energy F c of the system in terms of the one particle spectrum. 
The advantage of the canonical formulation is that statistical-mechanical entropy S c can 
be derived from F c by the standard rules. 

The second is the covariant Euclidean formulation. It starts with the calculation of 
the one-loop effective action W on the Gibbons- Hawking instanton. The free energy F E 
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in this approach is F E = where j3 is the inverse temperature. The covariant 

Euclidean formulation^ is especially useful in application to black hole thermodynamics. 

These two formulations are logically different, and in the presence of a black hole their 
comparison is non-trivial. The first part of the review is devoted to this problem. In 
Section 2 we define F E and F c on spacetimes without horizons and show the equivalence 
of these functionals. In Section 3 we discuss the features of quantum systems related to 
the horizon. Section 4 reviews the results of the canonical formulation for the spacetimes 
with horizons. Special attention is payed to divergences connected with the presence of the 
horizon and methods of their regularization. The covariant Euclidean formulation in the 
presence of a horizon is given in Section 5. We discuss ultraviolet one-loop divergences and 
show that the divergent parts of F E and F c are identical in the covariant regularizations, 
such as Pauli-Villars and dimensional ones. 

In the second part of the review we discuss the relation between statistical-mechanical 
entropy and the observable thermodynamical entropy of a black hole. In Sections 6,7 
we demonstrate how the divergencies of statistical-mechanical entropy are removed by 
the renormalization of gravitational couplings in the tree-level black hole entropy. The 
presence of the bare classical entropy make it impossible to give a statistical-mechanical 
explanation of the black hole entropy in the theories which require the ultraviolet renor- 
malization. This difficulty does not exist in the theories of induced gravity. The mech- 
anism of generation of the Bekenstein-Hawking entropy in such theories is described in 
Section 8. Our summary and conclusions are represented in Section 9. The Appendix 
is devoted to interpretation of the Noether charge which appears in the renormalization 
formula for the entropy because of nonminimal couplings of the fields. 

2 Statistical mechanics of quantum fields in a static 
gravitational field without horizons 

2.1 Description of the system 

We begin with discussion of statistical mechanics of quantum fields in a static gravitational 
field without horizons. This theory has been formulated and investigated in a number 
of papers starting with pioneering works by Gibbons pT]|,|]52]|, Gibbons and Perry JSSfl , 
Dowker and Kennedy |51| in the early seventies. A typical example is a quantum field at 



finite temperature in the gravitational field of a static non-rotating star. The gravitational 
field is described by the metric 

ds 2 = g 00 dt 2 + g ab dx a dx b , a, b = 1, 2, 3 , (2.1) 



2 The terms "canonical" and "covariant Euclidean" reflect the form of the presentation of the corre- 
sponding partition functions (in form of summing over the energy levels and as a functional integral, 
respectively). It should be emphasized that there is no standard terminology. We use superscripts C and 
E to refer to the quantities calculated in the canonical and covariant Euclidean formulations. 
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where goo(x) < 0. The metric ( |2.1| ) depends only on the spatial coordinates x a , and so it 
admits the Killing field ( = d/dt. On the spatial infinity the background is asymptotically 
flat and, by assumption, a time component of the metric goo tends to — 1. Because metric 
( |2.1D does not depend on time, one can easily define the statistical-mechanical ensem- 
bles of different fields on this background. We will be dealing with canonical ensembles 
characterized by the temperature T = measured at asymptotic infinity^. The local 
Tolman temperature measured by an observer at a point x a is T\ oc = | (700 1 1//2 /3 1 . 

To investigate the cases of Bose and Fermi statistics we will be considering, as an 
example, free scalar (0) and Dirac (ip) fields. The fields obey the Klein-Gordon and Dirac 
equations 

(- V^V M + m 2 + = , (2.2) 
( 7 ^V M + m)V = , (2.3) 

respectively. Here R is the scalar curvature and £ is the parameter of the non-minimal 
coupling. The covariant derivatives V M are defined according with the spin of the field. 
The Dirac 7-matrices 7 M = (7 , 7°) are defined by the standard relations { 7 M , 7^} = 2g fMU . 
Note that 7 is anti-Hermitean matrix. We define the spinor derivative as V M = d^ + T^, 
where = |[7 A ,7 P ] V^V^Va is the connection and Vj; are the tetrads. 



The following scalar products defined for solutions of equations ( |2.2j ) and (2.3) 



< 



fa, fa >=i J b d 3 x ^booh 1 {ftdtfa-fadtfi) , (2.4) 
<fa,fa>= [ A/%^ 2 , (2.5) 



where ^g = det g a b, are independent of the choice of the total Cauchy surface B. 

2.2 Canonical formulation and single-particle spectrum 

A canonical ensemble at temperature (3~ l is determined by the partition function 

Z c {(3) = Tr e-^ : . (2.6) 

The operator : 7i : is the Hamiltonian of the secondary quantized field. It determines a 
unitary evolution of the quantum field with respect to the Killing time t. As usual we use 
the decomposition of field operators onto positive and negative frequencies with respect to 



3 It is well known that an infinite bath of thermal radiation is a gravitationally unstable system. In 
order to deal with a stable situation one may consider a cavity of a finite size filled with the radiation. To 
characterize the system one can fix the temperature at the boundary or use the red-shifted temperature 
at infinity. These remarks are also valid for the case of a black hole being in thermal equilibrium with 
the radiation. For a non-rotating black hole of mass M the size of the cavity must be less than 3M S] . 
Speaking about the canonical ensemble we always assume that such a boundary exist. Our main subject 
is properties of ensembles of quantum fields connected with the presence of the horizon. The presence of 
external boundaries is not important for this consideration. For this reason we are not discussing them 
in more details. 
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the time t in order to define the creation and annihilation operators. The normal ordering 
in ( j2.6|) is with respect to these operators. For this normal ordering the energy of the state 
with zero temperature (vacuum) vanishes. If the vacuum energy is non-trivial, definition 
(|2.6|) can be easily modified. We will discuss this modification later. The canonical free 
energy of the system is 

F c ((3) = -(3- 1 \nZ c (f3) . (2.7) 

To proceed with the computation of F c {j3) it is convenient to rewrite Eq. (|2.7[) in 
another equivalent form based on a single-particle spectrum. Let u be frequency of a field 
mode with respect to the Killing time t. We call the set of these frequencies the single- 
particle spectrum. The spectrum of u is uniquely defined by the boundary conditions 
imposed on the system. If the system is given in a finite region with Dirichlet or other 
conditions on the boundary the spectrum is discrete. Some frequencies in this case can 
coincide and we introduce the corresponding degeneracy factor d(u). Then Eq. ( |2.7|) can 



be identically rewritten in the form [56] 



F c [(3] =r ] p- 1 Y,d(u;)\n(l - rje'^) . (2.8) 

The factor rj is related to the statistics, and it takes values t] = 1 and rj = — 1 for Bose 
and Fermi fields, respectively. Equation ( |2.8j ) is well-defined: the degeneracies d{uj) of 
the three-dimensional elliptic operators grow as u 2 at large lj but due to the exponential 
cutoff series ( [2.8]) converges. 

When the system has the infinite size the spectrum of u is fixed by the appropriate 
asymptotic conditions. A usual requirement is that fields fall fast enough at the spatial 
infinity. In this case the spectrum is continuous and the sum in ( |2.8| ) has to be replaced 
by the integral 

F c \J3]=r J l3- 1 (™ du d ^-\n(l-r l e-^) . (2.9) 
Jo aw 

The quantity dn ^ du is the number of levels in the interval (u,u + du). Equation (|2.9| ) 
can be obtained from ( j2.8| ) in the limit when the interval between the levels shrinks to 
zero. 

2.3 Single-particle Hamiltonian 

The single-particle spectrum can be found from the wave equations of the fields, Eqs. 



2D, ( |2.3| ). On the static space ( |2.1| ) these equations can be rewritten in the "3 + 1" form 



(dt + Hi)<P = , (id t -H d )^ = , (2.10) 



where the subscripts s and d are used for scalar fields and Dirac spinors, respectively. Hi 
and Hd are three-dimensional differential operators 

H 2 = \g 00 \(-V a V a -w a V a + m 2 + ^R) , (2.11) 
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7 a (V a + -w a ) +m 



(2.12) 



The 3-dimensional covariant derivative V a is defined in terms of the metric g ao on the 
hypersurface of constant time t = const. We denote the constant time surface as B. 
Operations with the index a are performed with the help of 3D-metric g a \, on B, see 
Eq. (|2.1|), and w a = |V a ln |g 00 | is the three-dimensional part of the acceleration vector 
= (0, w a ) of the Killing observer. 

The operators H s and Hd are quantum-mechanical single-particle Hamiltonians be- 
cause their eigen- values determine single particle spectra. By substituting the wave func- 
tions with the fixed energies <j>(t,x) = e~ laJ *0 a; (x), ip(t, x) = e _laJ V^( x ) into Eqs. ( |2.10| ) 
we arrive at the eigen- value problems 

#X,(x) = o; 2 ^(x) , Hfy u (x) = u 2 if> u (x) . (2.13) 

It is easy to check that H s and Hd are Hermitean operators with respect to the following 
inner products 



I d 3 x ^IsoohVKx) 2 (x) , (2.14) 

■IB v 



(Vi,^) = J b d 3 x ^^(x)^(x) , (2.15) 
where ?/> + denotes a Hermitean conjugated spinor. The form of relations (|2.14 ) and ( |2.15| ) 



follows from inner products ( |2.4| )-( |2~5| ) for four-dimensional fields. The inner products 
( |2.14| ) and ( |2.15| ) are used to normalize the modes. 

For the convenience of the computations we will also use another representation of H 2 
and Hd which can be obtained by the following transformation of functions and operators 

= e-> , ^ = e"5> , (2.16) 

H 2 s =e-°Hy , H d = e-I° H d e-*° , (2.17) 

where a = — | In |goo|- The inner products for transformed functions can be obtained from 
Eqs. ( [2.14|) and (|2.15| ) and they have the universal form 

($i,l> 2 )= / J®!j&x ($i) + $ 2 . (2.18) 

JB v 

Here $ denotes either scalar or spinor field, $>g = detg ao , and g ao = g a b/\9oo\ = ^dab- 
Let us stress that transformations (|2.16|) and Q2.17| ) do not change the spectra determined 
by Eqs. (|2.13|) . Thus, the operators Hi and Hi {i = s,d) are equivalent. 
From Eqs. (|2.17|) one finds the spinor Hamiltonian 

H d = *7o(7 a V a + e-°m) , (2.19) 

where {7^, 7^} = 2g^ u and g^ v = 57^ /| Pool- Thus, for the square of scalar and spinor 
operators we have 

H 2 = -V a V a + e- 2a m 2 + V { , i = s, d . (2.20) 
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The derivatives V a are defined in terms of metric g a b on the three-dimensional hypersurface 
B. The potential V{ is determined by the geometry of the background spaces and by the 
acceleration of the Killing observer 

V s = £R+ e~ 2CT (l - 60 (VX - uAu M ) , (2.21) 
V d = ^R + e~ 2(7 mYw^ , (2.22) 
R = e - 2<T [R + 6(VX - w^w^)] . (2.23) 

Note that the operators Hi can be found if transformation ( |2.16| ) of the fields is applied 
to wave equations (|2.2|) - (|2.3|) . The new equations obtained in this way can be expressed 
as equations of the fields on the ultrastatic spacetime with the metric 

ds 2 = -dt 2 + g ab dx a dx b , (2.24) 

related to the physical metric by the conformal transformation, g^ v = g^/l^oo!- The 
scalar curvature R for this metric is ( f2.23j ). Then developing the (3 + l)-formalism gives 



single-particle Hamiltonians (p. 20 ). In what follows all quantities calculated with respect 
to the ultrastatic metric g^ v will be denoted with a bar. The single-particle spectra and 
canonical formulations of the conformally related theories are equivalent. It should be 
emphasized that in the general case (in the presence of mass m and non-minimal coupling 
£ 7^ 1/6) the theory is not conformal invariant. For scalar fields the conformal invariance 
occurs only in the case when m = and £ = 1/6. 

2.4 Covariant Euclidean formulation 

The canonical formulation of statistical-mechanics is given in accordance with the unitary 
evolution of the system along the Killing time. It is explicitly related to the "3 + 1" decom- 
position and, therefore, it is not explicitly covariant. The covariant Euclidean approach to 
quantum fields at finite temperatures on stationary backgrounds was suggested by Gib- 
bons and Hawking pSfl , [j5"9] . This approach proved to be especially useful in application 
to thermodynamics of black holes ||58||-||61]|. 

Consider a manifold A4p with the Euclidean metric 

ds 2 = g TT dr 2 + g OL pdx a dx p , < r < (3 , (2.25) 

which is obtained from the static Lorentzian metric (|2.1| ) by the Wick rotation £ — > r = it, 
9tt = \goo\, and imposing the periodicity condition on the imaginary time r. We assume 
that the space is asymptotically flat and has the topology 1R 3 , so that the topology of 
Mp is H 3 x S 1 . 



According to Gibbons and Hawking [58], [59], the partition function Z and the effec- 



tive action W for a canonical ensemble of the fields $ in an external static gravitational 
background are defined by the path integral 

Z E ((3) = e ~ WM = J [D^e' 1 ^ . (2.26) 
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Here I[g, $] is a classical Euclidean action on the manifold M.p. Since g TT = 1 at spatial 
infinity, the parameter (3 is the length of S 1 and hence it has the meaning of the inverse 
temperature measured at the spatial infinity. Fields $ can have either Bose or Fermi 
statistics. Bose variables are assumed to be periodic in Euclidean time r with the period 
(3, while Fermi fields are antiperiodic. [-D$] is a covariant integration measure. For free 
scalar and Dirac fields the integration in ( |2.26|) gives 

W[g,(3] = W s [g,(3} + W d [g,P} , (2.27) 

W s [g,{3} = Uogd e t g - 2 L s , W d [g,{3} = -\ogdetQ- 1 L d . (2.28) 

The functionals Wi[g, f3\ are ultraviolet divergent and it is assumed that their divergencies 
are regularized, g is an arbitrary renormalization parameter with the dimension of the 
length, which does not depend on the background metric. In what follows we put q = 1 
for simplicity. If necessary the dependence of the effective actions on g can be easily 
restored. The operators Lj correspond to Eqs. (|2.2j ) and (|2~3| ) and read 



L s = -WV^ + £R + m 2 , L d = 75 ( 7 ^V M + m) . (2.29) 

Note that under the Wick rotation from Lorentzian to Euclidean signature the matrix 70 
has to be replaced by 270. The matrix 75 anticommutes with other 7's and is normalized as 
7g = 1. Both operators (|2.29|) are Hermitean with respect to the standard inner product 

($ 1 ,$ 2 ) = J d*xy/g$t®2 • (2-30) 

The Euclidean free energy F E [g,0\ is defined by the effective action of the system 

F l E [g,P} = p- 1 W l [g,P]-E^[g] . (2.31) 

Similarly to Ff[g,j3], the so-defined Euclidean free-energy Ff[g,(3] vanishes at the zero 
temperature. This makes the comparison of the two free energies more simple. The 
vacuum energy 

E*[g} = lim (p^Wfafl) (2.32) 

does not contribute to the entropy. Since the free energy F t E [g, (3\ and the vacuum energy 
Ef[g] are defined in terms of the covariant Euclidean action Wi[g, (3\ this approach is called 
covariant Euclidean formulation. 

2.5 Relation between canonical and covariant Euclidean formu- 
lations 

The relation between canonical and Euclidean formulations in the absence of the horizon 



was discussed by a number of authors |pT|-|M, |57|] , ||62|-||65|. 
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For the comparison of these two formulations the crucial role is played by the repre- 
sentation of the canonical free energy in terms of the effective action in the ultrastatic 
space M.p with the metric 



ds 2 = dr 2 + g a adx a dx p , < r < (3 



[2M) 



which is conformally related to M.p. This space is the product S 1 x B. 

Let us define the operators Li on M.p which are conformally related to operators Lj 

Eq. 



-3<T 



L, e a 



L d = e-^L d e* a 



It is easy to show that 



is = m 



di 



L d = w T (H d + d T ) , L 2 = H 2 -d 2 T . 
For these operators one can define the effective actions 

W s [g, 0\ = - log det L s , W d [g,(3] = - log det L d 



(2.34) 

(2.35) 
(2.36) 

(2.37) 



(For convenience we consider W d as the functional of the physical metric g^.) 

For the system with a discrete spectrum one can find with the help of Eqs. fl2.35p and 
( [2.36 ) the following relation between the canonical free energy and the effective action in 
the ultrastatic space is |57| 

FC{g,f3}=f3- 1 W l {g,P}-E°ig} . (2.38) 
The quantity E^[g] is the vacuum energy for the fields in the ultrastatic space 



to 



^2.39) 



Eqs. ( |2.38| ) and Q2.39| ) can be generalized to include the systems with continuous spectra. 
For this purpose one has to take the limit when the intervals between the levels vanish 
and to replace sums over u with integrals. The derivation of ( |2.38| ) and (|2.39| ) can be 



found in work by Allen [0], see also |6( 

It is important to note that Ef[g] includes all ultraviolet divergences of the functionals 
Wi[g, (3\. Thus the functional Ff[g, 0\ is ultraviolet finite. The geometrical structure of the 
divergences does not depend on the temperature of the system j54j. This is a consequence 
of a more general property of the quantum field theory whose ultraviolet singularities do 
not depend on the quantum state of the system f67| . The renormalization of Wj[<7,/3] is 
equivalent to renormalization of 

Equations ( 2.34 ) are crucial for finding relation between canonical and Euclidean free 
energies. The classical actions corresponding to the operators Lj and Lj are 



I i\9^i\= d 4 x ^/g ip+ L^i , lf[g,cpi)=l d A x Jg (pf 

JMp JM/3 



(2.40) 
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where the notation (ft is used for scalars or spinors if). So, as a result of Q2.34| ), 



(2.41) 



for (f) = e~ a (J) an d ^ = e~^ a ip. In case of massless scalars with £ = | or massless spinors 
fields the operators Li and Li have the same form which means that classical theories are 
conformally invariant. In general case this invariance is broken. However, as it was shown 
by Dowker and Schofield [ 64 1 , |p5|1 it is still possible to introduce an auxiliary conformal 



charge in the classical actions and interpret Eq.( [2.4l|) in terms of a pseudo conformal 
invariance. 

In quantum theory the classical symmetries are known to be broken because of anoma- 
lies []67l,|l68|. Thus, the relation between the renormalized actions has the form 



Wi[g,P} = Wi[g,(3] + (3ni[g} 



[2A2) 



For scalar and spinors fields the anomalous terms were found explicitly in p3,[BS 

The anomaly is proportional to (3 and so it contributes to the vacuum energy only 



E^la] = ma] + ^[g] 

As a result, the free energies F® and Ff, coincide 

Ff\g,0\ = Ff\g,0\ 



[2A3) 



[2.U) 



Therefore, on the static spacetimes without horizons the covariant Euclidean and canon- 
ical formulations of statistical mechanics of quantum fields are equivalent. 

The appearance of the anomalous terms in ( |2.42| ) can be attributed to non invariance 
of the integration measure with respect to conformal transformations |j!J-[[7[J. Path 
integral definition ( |2.26| ) of functionals Wi[g, (5\ employs the covariant measures g x ^d(\> for 
scalars and g~ Vd ^ 2 dip + dip for spinors, where g = detg^ u , and is the dimensionality of 
the spinor representation (see for instance Refs.[[TT[]-[[72]). The difference in the prefactors 
in scalar and spinor measures is related to the different integration rules for Bose and 
Fermi variables. The integration measure for actions W^g,/?] is different. It is g^^d^ 
and g~ Td / 2 dip + d-ip, for scalars and spinors, respectively, and it is covariant with respect 
to the ultrastatic background. One can obtain this form of the measure by canonically 
quantizing the theory [f73[ . 



3 Features related to the horizon 

Many aspects of the computation of the statistical-mechanical entropy of a thermal at- 
mosphere of a black hole are directly related to the presence of the horizon. Consider a 
general static spacetime with the Killing vector ( = d t which has a Killing horizon where 
( 2 = 0. At the moment we do not require that the metric obeys the Einstein equations. 
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If the surface gravity 

1/2 

(3-1) 



-~(C^)l c *=o 



does not vanish, the metric near the horizon can be represented in the form 

ds 2 ~ -K 2 p 2 dt 2 + dp 2 + dfl 2 . (3.2) 

Here dVt 2 is the metric on the two-dimensional bifurcation surface of the horizons where 
Q L = 0. We denote this surface E. In coordinates ( |3.2|) the horizon is located at p = 0. 
Such form of the metric is general for non-extremal black holes. 

In the presence of the horizon the spectrum of quantum fields with respect to the 
Killing time has a number of important new properties: 

1. The single-particle spectrum is continuous even if the system has a finite size, i.e., 
when there are boundary conditions imposed on the system at the finite distance 
from the horizon; 

2. Regardless of the spin and mass of the fields, the spectrum of uj runs down to uo = 
and, thus, the usual mass gap is absent; 

3. The bifurcation surface S is invariant under the time-evolution. 



The continuity of the spectrum can be easily understood in representation ( 2.20 ) for 



the single-particle Hamiltonians. The operators Hi are given on the spatial part B of 



the ultrastatic space ( |2.24|) . The space B is always non-compact because the conformal 
transformation is singular at the event horizon and it moves the points of £ to the spatial 
infinity in B. 



The mass gap for the operators Hi vanishes for the following reason. Equation (|2.2C| ) 
shows that the mass rrii of the field does not have any effect near the horizon because 
of the factor e~ 2a = \goo\- The operator — V a V a has a mass gap which is of the pure 
geometrical origin. It arises since the space B has asymptotic constant curvature near 
the horizon. However, on the horizon V s = —k 2 , Va, = — § n 2 . These potentials act as a 



tachionic mass |74|,|[75| which exactly cancels the mass gap caused by the curvature of B. 

Thus, the system we are dealing with behaves near the horizon similar to a massless 
quantum theory on a non-compact space. It is well known that such theories run into the 
difficulties related to the infrared divergences. Because of these divergences the densities 
of eigen-values of quantum-mechanical Hamiltonians Hi blow up near the horizon. As a 



result, free energy (|2.9| ) in the canonical formulation diverges at any temperature. 

In the presence of the horizon the covariant Euclidean formulation also exhibits new 



features. In this case the Euclidean section Mp, Eq. ( |2.25| ), of the black hole solution 
cannot be regular for arbitrary values of /3. Near the surface £ where the Killing vector 
( = d T vanishes (the Euclidean horizon) the metric has the form 

ds 2 ~ n 2 p 2 dr 2 + dp 2 + dn 2 , < r < (3 , (3.3) 
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which follows from Eq. ( |3.2| ). Thus, near this surface Aip looks as Cp x S, where Cp 
is a conical space. It is easy to see that the conical singularity disappears and Aip is 
regular space only at the special value P = Ph = 27r/t~ 1 . The corresponding temperature 
T H = Pjj 1 coincides with the Hawking temperature and the corresponding quantum state 
is known as the Hartle-Hawking vacuum J76| . The physical meaning of Th is that it is the 



temperature of the Hawking quanta emitted by an evaporating black hole. T# also gives 
the temperature at which the quantum radiation can be in thermal equilibrium with a 
black hole. 

For P ^ Ph conical singularities result in additional ultraviolet divergences in the ef- 
fective action on Aip. Thus, in the presence of the horizon both canonical and covariant 
Euclidean formulations acquire new divergences. These divergences are of the different ori- 
gin; they are infrared in the canonical method, and ultraviolet in the covariant Euclidean 
approach. In Sections 5 and 6 we develop canonical and covariant Euclidean formulation 
in the presence of the horizon and establish the relation between the divergences in these 
formulations. 

4 Canonical formulation in the presence of a horizon 
4.1 Density of levels and its properties 

The divergence of statistical-mechanical quantities in the presence of a horizon is directly 
related to the infinite growth of the density of states of the Hamiltonians. To investigate 
this property we follow Ref . |J?J . The idea, which is close to the earlier approach of , 
is to relate drii(uj) / du to the heat kernel of the operator H 2 . The latter is an elliptic 
operator and its heat kernel is well known. For a continuous spectrum one has 

Tre-^= r d J^e-^ . (4.1) 
Jo duj 



The density drii(uj)/duj can be found from (|4.1| ) with the help of the inverse Laplace 



arc 



transform [fT7|| . The diagonal matrix elements (x\exp(—Hft)\x) = exp(-Hft)^^ 
well defined and finite. However, the corresponding trace which involves the integration 
over the noncompact space & diverges. To understand why it happens it is sufficient 
to study the behavior of exp(-Hft) ^ near the horizon. To estimate the leading 
asymptotics we can neglect the curvature of the two dimensional surface of the horizon 
and approximate the black hole metric (|3.2|) by the metric on the Rindler space 



ds 2 = —n 2 p 2 dt 2 + dp 2 + dz\ + dz\ , — oo < zi, z-i < oo , p > . (4.2) 
Then the metric on the conformal space B 

dl 2 = K~ 2 p~ 2 (dp 2 + dz\ + dz\) (4.3) 
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coincides with the metric of the hyperbolic manifold IH 3 of constant negative curvature 
R = —6k 2 . A review of the Laplace operators and their heat kernels on hyperbolic spaces 
can be found in Refs. J78[ and |79| . The eigen-functions of Hf are completely determined 
by the requirement to have the correct decay properties at infinity p —>■ oo |f75| . Thus, 
no additional conditions at the horizon p = are needed. When p — > the fields become 
effectively massless and the diagonal elements of the scalar and spinor heat kernels on Hi 3 



are known exactly |7 
-tm 



diag (4vrt) 3 / 2 



-tm 



diag (4vrt) 3 / 2 



1 



-kH 



(4.4) 



A summation over the spinor indexes is assumed and it gives the factor r^. The geometry 
essentially differs from the Rindler one far from the horizon and the mass term becomes 
important. For this reason in general there are corrections to Eqs. ([4.4|) proportional to 
the powers of p 2 . As it was shown in [66], the structure of these terms can be analysed 
by using asymptotic properties of the heat kernels on B. 

Therefore, as follows from Eqs. (|4.1|) and ( |4.4|) , the trace of the operators and their 



density of states grows as the volume of B. 

There are several ways how to regularize these divergences. For instance, one can 
restrict the spatial size of the system. In the given case it means that the region of the 
physical spacetime where the proper distance to the horizon is smaller than some length, 



say e, has to be excluded from the consideration. To this aim t'Hooft [0 suggested to 
impose the Dirichlet boundary conditions on the fields at a surface located outside the 
horizon and at the proper distance e from it. t'Hooft's approach is known as the "brick 
wall" model. A similar but simpler procedure, called the volume cutoff method, was 
proposed by Frolov and Novikov [[31]] (see also Refs. [p0[],||8l|). In this method all spatial 
integrations are cut off at the proper distance e without imposing a boundary condition. 

In the volume cutoff method one effectively cuts a region near the horizon which 
makes the spacetime incomplete. There exist other types of regularizations which allow 
one to work on the complete spacetime background. An example is the Pauli-Villars (PV) 
regularization which was first used for the problem by Demers, Lafrance and Myers [ j39fl . 
The density dn/du turns out to be finite even in the limit e — > but it depends on the 
regulator mass /i and for p, — > oo it grows as p? . Another option suggested in |48| is to use 
the dimensional regularization. The idea of this method is that the power of the leading 
divergency in Eqs. (|4.5| ) and fl4.6| ) depends on the number of spacetime dimensions. In 
/^-dimensional spacetime the leading divergence of the volume of B is e 2 ~ D , if D ^ 2. One 
can use Dasa regularization parameter for density of states, and take the limit e — > at 
Re.D < 2. The density dn/duj then has a pole at D = 4. 

In the volume cutoff method the divergence of dn/duj is infrared. In dimensional 
and PV regularizations the divergence of dnjdio can be directly connected with standard 
ultraviolet divergences. For this reason, we can speak about infrared and ultraviolet limits 
of the theory depending on which type of regularization is used. 
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4.2 Infrared limit and volume cutoff 



Let us denote the regularized density in the volume cutoff method as drii{uj\e) / dw and 
investigate its asymptotics in the limit e — > 0. After integrating the heat kernels over the 
region p > e and using the inverse Laplace transform in ( |4.1| ) one obtains the regularized 
expression for the divergent part of densities of states p6| 



dn s {uj e) 




1 ({ 


doj 




47r 2 K 3 Jy, [ 




div 








dn d (uj\e) 






dui 

L J c 



CO 



1 l^i 
— - -V In - 



k 2 e 2 



(4.5) 




(4.6) 

The notation J s assumes that the integration is performed over the bifurcation surface of 
the horizons E, so that J E 1 = A, where A is the area of E. We imposed an additional 
cutoff / at large distance p. The quantity r d is the dimension of the spinor representation, 
so that r<£ = 4 for 4 dimensional Dirac spinors. Let (i = 1, 2) be two unit mutually 
orthogonal vectors normal to E, then P^ v = Y^=i is a projector onto a two dimen- 
sional surface orthogonal to E. The quantities V and Q which enter ( fL5| ) and ( |4.6| ) are 
defined as 

V = 21l-Q , Q = P^R^ , Tl = P lxv P Xp R pXvp . (4.7) 

Note that the leading divergence e~ 2 in Eqs. (|4.5| )- (|4.6| ) is already present in the 
Rindler approximation ( ^4.4|) . The mass and non-zero curvature result in the additional 
logarithmic divergences ln(e 2 // 2 ). 

Substituting and Q into expression (|2.9|) for the canonical free energy we 

obtain 



2 

Qp^ 2 



7T 



720/3 4 A8(3 2 
7tt 2 



4' 



7vr 2 




+ 



U440/3 4 192/3 2 / e 2 



96/3 2 V8 12 



m 




(4.8) 



(4.9) 



5760/3 4 ' 

This enables one to calculate other characteristics of canonical ensembles. In particular, 
one can find the e divergence of the statistical-mechanical entropy 



t dF°Jg,/3,e] 



d/3 



(4.10) 
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Evaluated at the Hawking temperature (3 1 = P H X = k/2it the divergent part of the 
contributions of bosons and fermions to the entropy is 



f_L_ 

/s \360e 2 



1440 



2TZ - Q + 30 



Hi? - 30m 2 



In 



/ 2 



5 fdiv[^.^H,e] 




1 



1440e 2 5760 



7n + AQ-5R-Q0m 2 




(4.11) 



(4.12) 



As we already mentioned, the vacuum energy omitted in (|2.6| ) does not contribute to the 
entropy. For scalars, the leading divergence e~ 2 of the entropy is determined by j3~ 4 term 
in free energy (f4.8|). For the spinor fields the situation is different. In order to find e~ 2 



divergence to the entropy one has to know both the leading, j3 4 , and subleading, f3~ 



terms in free energy ( f4.9|) . 



To determine the divergences of the free energy and entropy one has to know only 
the asymptotics of the heat kernel operators j6"6|| . This is not sufficient if one wants to 
calculate not only the divergencies but the quantities themselves. Important examples, 
when dn/du, F c and S can be calculated exactly, were studied in Refs. 



The authors considered the case when the spatial part of the spacetime is of the form 
B = IR + x E, where X is a manifold of constant curvature. In particular, for the Rindler 
spacetime S is a two dimensional plane and the density of states of massive scalar fields 



in Rindler space is [71 
dn s (u\e) A 



du 



4tt 2 k 3 



u 



+ 



m 2 K 2 . e 2 m 2 
In 



m 2 K 2 



1 + 2Re if>(iu/K)) 



(4.13) 



Here ip is the logarithmic derivative of the T-function. The terms which vanish as e — > 
are omitted. It follows from (|4 .13 ) that the role of additional cutoff I in ( |4.5|) and ( |4.6|) is 
played by wT 1 . 



4.3 Ultraviolet limit 

A. Dimensional regularization 



The dimensional regularization is the simplest scheme which enables one to define dn/du 
on the complete spacetime. In this regularization one can put e = from the very 
beginning. The quantity dn/du depends on the complex parameter D associated with 
the dimensionality of the spacetime and it has a pole singularity when D = 4 IBT 



dn s (u\D) 
du 



dn d (u\D) 
du 



D-4 



r(i-f 

(4tt)^/ 2 k 



m 2 — 



r(i-f 

( 4vr )D/2 



III 



D-4 



2 M R Q 

m H 

12 8 



V 

K 2 



— V 



(4.14) 



(4.15) 
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Table V. 



spin 


di 


h 


Qi 


<?2 


93 


Pi 


P2 


P3 





J"* 


1 


1 


-1 


|(1 -60 2 


l 

60 







l 

2 






-|r d 




td 




k r d 


~TE r d 



These expressions can be used to get the divergences of the canonical free energy for 
scalars and spinors 



Fg[g,P,D] = - V 



r i 



D 



2 TT "m 



2^D-4 



(4tt) d / 2 3k(3 2 
The corresponding divergence of the entropy is 

S^ iv [g,P,D} = p 2 

r(l-f) 27r 2 m D - 4 



/ 4tt 2 

fim 2 - Pi-^V + p 2 R + p 3 Q) 



K 2 (3 



(4.16) 



dF^\g,P,D] 
dp 



(4tt) d / 2 3nf3 



57T 



A m ~ I Pi T^2 P +P2R + P3Q) 



(4.17) 

Constants /1 and depend on spin and are given in Table 1. This table contains also 
other similar constants which we meet later. 



B. Pauli-Villars regularization 

It is known that the dimensional regularization reproduces only the logarithmic diver- 
gences. For this reason it is also worth studying F^ iv in another, more complete regular- 
ization. For our purpose it is convenient to use the Pauli- Villars method. In this method 
for each of the physical field, one introduces 5 additional auxiliary fields: 2 fields with 
masses which have the same statistics as the original field and 3 fields with masses M' r 
which have the wrong statistics^. To eliminate the divergences the masses of the auxiliary 
fields must obey the two restrictions 

/(l)=/(2) = , (4.18) 

where 

f{p)=m 2 P + J2M 2 k p -Y,(M^ p = . (4.19) 

k r 

4 They are fermions for scalars and bosons for spinors. 
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These constraints can be resolved by taking M 12 = y/3fi 2 + m 2 , M[ 2 = y/ fj, 2 + m 2 , Mg 
^4/i 2 + m 2 (see [j39| ). The regularized density of states in this method is 



drii(uj\fjL) _ dni(u,m) 



du 



du 



E 



dnj(uj,M k ) 
du 



E 



d ni {u,M' r ] 
du 



(4.20) 



The quantities drii(u , M^) / du and drii(u , M' r ) j 'duj are the density of states of the Pauli- 
Villars partners. The fields with the wrong statistics give negative contribution in the 
regularized density. Because the number of such fields equals the number of the fields with 
the proper statistics the leading e-divergences in relations (|4.5|) and (|4.6|) are cancelled. 
Logarithmic divergences lne 2 are also canceled because of the constraint f{jp = 1) = 
imposed on masses. As a result, the regularized density of states ( |4.20|) does not contain 
the divergences when e — > and it can be defined on the complete background. Since in 
the presence of regulators dni(u\/j,) j 'du is finite one obtains the same answer as when one 
uses the dimensional regularization for its calculation. Using Eqs. ( 4.14 ) and ( 4.15Q and 
by taking into account constraints Q4.18|) one gets 



dn s (uj\fi) 
du 

dn d (u\n) 
du 



(4tt) 2 k 



26 + a 



u 



K 



V + 2 



l 



(4vt) 2 k 



R Q s 
2b + a\'—V-- + 4 



u 



K 



6 4 



The functions a and b depend on m and [i 

dp 



p=0 



In m 2 - ln M k + E HK) 

k r 



b = 



df 



dp 



r/\2 



(4.21) 
(4.22) 

(4.23) 
(4.24) 



= m 2 ln m 2 + £ M 2 ln M 2 - £(M;) 2 ln(M; 

p=l k r 

In the Pauli-Villars method /i 2 plays the role of the ultraviolet cutoff. In the limit of 
infinite masses M' r , i.e. at // — > oo, dni{iu\\i) j du is ultraviolet divergent. In this limit 



ln 



\i ln 



729 



2, V 

mm — : 



(4.25) 



m* 256 m A 

Thus, in general dn/du includes both quadratic and logarithmic divergences. 

From Eqs. (|4.21|) -( [4.22|) one can derive the divergence of the statistical-mechanical 
free energy 



V 



f 4tt 2 

6 /i + a [ Pl ^Jp V + P2R + P3Q 



(4.26) 



A special case of expression (|4.26|) for scalar fields on the Reissner-Nordstrom black hole 
background was first derived in [35]. The divergent part of the entropy obtained from 
Q4.26D at the Hawking temperature is 

V 



SZ[g, fa, 14 = ^^ Wi + a (2piV + p 2 R + PsQ)} 



(4.27) 
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By taking into account Eq. ( f4.25|) we can rewrite this expression at large fi in the form 



r 2 

487T je v ' m z 



(4.28) 



where c = ln|||. Expression ( |4.28| ) has the same structure as the entropy divergences 
( [4.11 ) and ( [4.12 ) in the volume cutoff regularization. It is easy to see that in these 
regularizations parameters e and /i -1 correspond to each other. The fact that Pauli-Villars 
regularization results in a cutoff of the integrals near the horizon at the proper distance 
comparable to the inverse masses of the fields allows the following interpretation. Near 
the horizon where the local temperature becomes greater than /z, the massive regulators 
are thermally excited and because of constraints ( [4.1 8|) their contribution exactly cancels 
the contribution of the physical field. 

It is interesting that for each field one can find a relation between e and by 
equating the leading divergencies. After that it is possible to find the connection between 
I with m _1 which makes equal logarithmic divergencies as well. Note, however, that in 
this identification the relation between e and [i is different for fields of the different spins, 
see Eqs. ( pp , fl4l2D and ( jgg ). 



4.4 WKB computations and the "brick wall" model 

The divergences in statistical-mechanical quantities can be also obtained by using WKB 
method. This way of computations was suggested by t'Hooft [PS] and then used by many 
authors, see for instance, Refs. |8"2T|-|]8~7|]. To illustrate this method we consider a scalar 



field on asymptotically flat spherically symmetric black-hole background 
d s 2 = -g(r)dt 2 + g-^^dr 2 + r 2 (d6 2 + sin 2 9 dip 2 ) . 



(4.29) 



Here r > r^, and is the horizon radius where g(rh) = 0. Schwarzschild and Reissner- 
Nordstrom black holes are described by metrics of this type. We are interested in modes 
( x ) whh energy u) and angular momentum i. They are solutions of the eigen- value 
problem 

H 2 ^^^ 2 ^ . (4.30) 
This equation is reduced to one-dimensional problem 

,2 



g 



»-dAr 2 d r )+V* f {r)-uj 



with the potential 



Veff(r) 







^g>) 2 + g (r)[m 2 + r~ 2 £(£ + l) 



(4.31) 



(4.32) 



Near the horizon this potential is negative K//( r ^) 
and Veff(j' = 00) = m 2 . 



-k 2 , while at r ^> it is positive 
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Let us consider now the "brick wall" model We assume that the "brick wall" 

Dirichlet condition is imposed at the proper distance e near the horizon, and r(e) is the 
location of the "brick wall" in coordinates ( |4.29| ). Equation (|4.31| ) enables one to estimate 



n s (uj\e 

TT 



in the quasiclassical approximation the number of energy levels n s (uj\e) with the energy 
smaller than uj 

Here £(u) is the maximal angular momentum at which the square root in ( J4.33j ) vanishes, 
Tg 3> Th is an additional infrared cutoff. 

To estimate the asymptotic behavior of n s (uj\e) at e — ► 0, we assume that uj 2 is large 
compared to the curvature of the background. Then only the contribution of the large 
momenta is important and the sum over i can be replaced by the integral, which can be 
easily calculated 

/ I * 2 f r B r 2 dr 2 , 3/2 A f uj 3 3 2 2 , e 2 \ , AriA . 



In the last equality we put A = 47rr^, dp = g x l 2 dr and use for the metric (|4.29|) near the 



horizon the Rindler approximation ( fL2|) . It is easy to see that Eq. ( |4.34f ) gives the same 
density of energy levels as expression (|4.5|) for a scalar field with V = R = 0. 

Thus, at least for the leading divergences there is an agreement between the volume 
cutoff regularization and the WKB "brick wall" model. This agreement remains if one 
chooses the Newman boundary condition instead of the Dirichlet one. The numbers 



n a {u\e) for these conditions differ only by a numerical constant |33|. It should be noted 
that for the physical field the horizon is not a real boundary. From the mathematical point 
of view it means that the wave functions in Eq.( [4.30D with uj < m can be completely fixed 
by the requirement to decay fast enough at r > r/j and no boundary conditions at the 
cutoff length e are required. For this reason the volume cutoff method seems to be more 
appropriate than the "brick wall" approach. 

The WKB method can be also used to study the ultraviolet divergences. It was 
employed, for instance, in 3|| in PV regularization. Analogous computations were done 



in Refs. [B3fl and ||S5|| . Almost all WKB computations are restricted to the case of the 



scalar field (see, however, Refs. |87j and |88 



5 Covariant Euclidean formulation 

5.1 Volume cutoff and high-temperature expansion 

In covariant Euclidean formulation of the theory, the problems related to the divergences 
on the horizon can be treated in the same way as in the canonical formulation; namely, by 
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introducing either the volume cutoff or ultraviolet regularizations. Obviously, comparing 
the Euclidean and canonical formulations has sense only for the equivalent regularizations. 



In the volume cutoff regularization, one effectively cuts a region near the horizon. This 
makes the spacetime incomplete. As the result, one gets a theory on a static spacetime 
without the horizon. As we discussed earlier, the covariant Euclidean and canonical 
formulations are equivalent in this case. 

The form of the Euclidean action near the horizon can be found by using the high 
temperature asymptotics. These asymptotics give a good approximation because the 
local temperature infinitely grows when approaching the horizon. The high temperature 
expansion of an effective action was obtained by Dowker and collaborators in [54], p4| 
and p5| (see also Ref. f89|). The expansion has the form 



(5.1) 



Wi\g,l3) = /3 f d 3 xg 1/2 [bfap) + h^x)] + w{%] - AWi 



For spinors Wj [#] = 0, while for scalars Wj^lg] = — ^('(0\H 2 ) is determined by the 
^-function of the operator H 2 , see |f)4|,|l65"||. The quantities 6, and hi which enter relation 
ni) are 



b s (x,(3) 



7T" 



1 



90$ 4 24/52 



m 



+ 



720/3? 48/3) 



-^R + ^(Vu; - w 2 ) + m 2 



a s ,2(x) lift 
167T 2 2n 



+ 



, In- 

16tt 2 2tt 



h s - 



1 



2880vr 2 
r r j 



5w 2 (w 2 - 2Vw) - 3(Vw) 2 + R^wW - 30 [(£ )R + m 
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w 



7w 2 (w 2 - 2Vw) - 18(Vw) 2 + TlR ilv w il w v - 5(R + 6m 2 



1440tt 2 

In these expressions w 2 = w^w^, Vw = V^Wu, @i(x) 



(5.2) 
(5.3) 



(5.4) 
(5.5) 



i7oo 1 1//2 /3 is the inverse local Tolman 



temperature, a^x) are the second heat coefficients of the 4-dimensional operators Li 
Eq.( |2.29| ). The trace over spinor indexes in a d)2 (x) is assumed. 
The quantity AWi which enters (|5.1|) is 



AWi 



167r 5/2 



n=3 



Ci,nT (n - - ) C, R {2n - 3) a 




2n-4 



Here c s>n = 1, Cd, n = 1— 2 2n 3 , (r(z) is the Riemann ^-function f90fl . are the coefficients 



of asymptotic expansions of the heat kernels of the operators H 2 . The quantity AWi can 
also be rewritten in terms of the physical metric and local temperature [|64 l,[l65|. Let us 
emphasize that except the term W^ the high temperature expansion ( |5.1| ) of the effective 
action has a local form. The non-local contributions to the finite-temperature effective 



action was studied in |H[ . 
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In the volume cutoff method the integration in ( |5.1| ) is performed over the spacetime 
of the black hole till the proper distance e from the horizon. The action Wi[g, f3,e] and 
free energy Ff[g, (3, e] diverge when e — > 0. The divergences are caused only by the terms 
proportional to /3~ A and (3~ 2 in the functions bi(x,/3). The functions hi contribute to the 
divergence of the vacuum energy. As expected, the divergences in Euclidean and canonical 
free energies coincide 

F t ? div [g,P,e] = F* w [g,P,e] . (5.6) 

Here the quantities Ff div \g, (3, e] are determined by Eqs. (|4.8|) , ([4~9D and Ff dlv [g, (3, e] can 
be found from the actions (|Q|) by subtracting the vacuum energy, see Eq. ( [2.31| ). 

5.2 Conical singularities and ultraviolet divergences 

We now consider the definition of the covariant Euclidean effective action on a complete 
manifold Mp, i.e., on a manifold with conical singularities^. The important feature of this 
case is that the conical singularities result in additional ultraviolet divergences depending 
on j3. 

It is convenient to introduce the wave operators Aj = — V M V M + X{, acting on the 
scalars and spinors. Here V's are corresponding covariant derivatives. We have 

L s = A s + m 2 , L 2 d = A d + m 2 , (5.7) 

where X s = (1/6 — £)R and X d = \RI. The analysis of Laplacian operators on cones and 
their heat kernels can be found in Refs. [[L00|l - |102|| . For physical applications the suitable 
representation for the Green functions and heat kernels of integer and half-odd-integer 
spins was given by Dowker [|103| |,| |104[ , see also Ref. ||105|| . These results are based on the 



generalization of the representation obtained by Sommerfeld 100 years ago [|106|1 . 

In general, a one-loop effective action can be defined by the Schwinger-DeWitt repre- 
sentation 

V f°° ds __ m W „- s a 



W[g,{3] = ^logdet(A + m 2 ) = -j- / — e~ m s Tr e" sA , (5.8) 

where 5 2 is an ultraviolet cutoff. (For briefness we omit the index i indicating what type 
of the fields we are considering.) The geometrical structure of the ultraviolet divergences 
is determined by the first terms in the asymptotic expansion of the heat kernel at small 
s. On manifolds without boundaries it takes the form 

Tre-^^^-^(B + sB 1 + s 2 B 2 + ...) , (5.9) 

where D is the dimension of The heat (Hadamard-Minackshisundaram-DeWitt- 

Seeley) coefficients B k for k > 1 can be represented as the sum of two terms 

B k = A k + A Ptk . (5.10) 

5 There are several explicit calculations of an effective action on some two-dimensional |f)2f , fsof . |9^] . 
three-dimensional [p4| and four-dimensional [p5| , |p6|] spaces with conical singularities. It is worth men- 
tioning that conical singularities appear in a large number of other physical applications including cosmic 
strings |97j , orbifolds in string theory |p8| and topological defects in superfluid phases of Helium 3 (99) . 
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Here has the form of the standard coefficient defined on the regular domain of J\Ap. 
The term Apk is an addition due to conical singularities. This contribution is a functional 
on E which depends on the geometrical characteristics of M.p near this surface. The first 
two coefficients A^ and Apj. have the form 



A l = d 1 R , Ap \ = —MY - l)A 

M = -^l UR^OR^xp + qiR^R^v + q 3 R' 
180 JM B -T, v 



A M = — 



l 0- 

37 JS 



( 7 4 -l)p 1 P + (7 2 -l)(p 2 i? + p 3 Q) 



(5.11) 

(5.12) 
(5.13) 



where 7 
Table 1. 



The coefficients di, fi, qi and Pi for spins s = and 1/2 can be found in 



For a scalar field the coefficient Ap^ was first found by Cheeger ||100|| , see also Refs 



107|| , ||108]1 . The coefficients Ap } % for spins 1/2 and 1 were obtained in Refs. [|109|| , ||TT0 



and for spins 3/2 and 2 in [ |110|| . The coefficient Ap^ and general structure of the higher 



coefficients Ap t k were analysed in Refs. [Ill] , ||112|| and ||113|| for scalars, in Ref. [B6|, and 
for spinors in Ref. [|TT^ 



The divergent part W^vl^f, /3] of the effective action on J\4p can be written in different 
regularizations. It is convenient to begin with the dimensional regularization. For 
one finds from Eqs. (|5.8|) and 



W div [g,P,D] = -^- [ 
I Jo 



00 ds 



s 2 B, 



(5.14) 



(4tts) d / 2 

The Euclidean free energy is obtained from the action by subtracting the vacuum energy 
(see Eq. ( |2.31|) ). On a regular Euclidean manifold the divergences are determined only 
by the coefficients A , Ai, A 2 and they are proportional to the period j3. So they do 
not contribute to the free energy and the entropy. In case of conical singularities the 
divergences have the form of polynomials in powers of (3~ l because of the additional terms 
Ap t /,. These terms are not removed by subtracting the vacuum part. For this reason the 
Euclidean free energy is divergent. The divergent part F^ v [g, (3, D] can be easily found 
with the help of Eqs. ( pTT| ), (pp and |U4|). 

As follows from Eq. (|4.16|) , the ultraviolet divergences of the canonical and Euclidean 
free energies for scalar and spinor fields are identical in the dimensional regularization. 
This coincidence also takes place in the Pauli-Villars regularization, so one can write 



s, d 



(5.15) 



where 6 is a regularization parameter (5 = D — 4 for the dimensional regularization and 
5 = yT x for PV regularization). As the important consequence of ( |5.15p , the divergences 
of the entropy in different formulations coincide as well 



s, d 



(5.16) 
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5.3 Relation between canonical and Euclidean formulations in 
the presence of horizons 

Till now we discussed and compared the divergent parts of the free energies calculated 
in the canonical and covariant Euclidean formulations in the presence of a horizon. In 
this Section we make remarks concerning the relation between the finite parts of the free 
energies in these formulations. 

As we discussed earlier, the canonical and Euclidean formulations are completely 
equivalent for static spacetimes without horizons. This conclusion is based on the fact that 
the effective actions W[<7,/3] and W\g, 0\ given on M.p and M.p, respectively, are related 
by a conformal transformation. The Euclidean and canonical free energies obtained from 
W[<7, /3] and VV"[<7,/3] by subtracting zero temperature parts do not have the divergences. 
They are free from the anomalies caused by renormalization and as the result F E and F 
coincide in the absence of a horizon. 

The total functionals F c [g,/3,e] and F E [g,/3,e\ (containing both e-divergent and reg- 
ular parts) also coincide on the backgrounds with horizons provided one uses the volume 
cutoff method. In this method the conformal transformation which relates the free ener- 
gies is well-defined. 

The difficulties arise when one compares two formulations on a complete background 
which includes the horizon. In this case the spaces Aip and M.p, have different topolo- 
gies, IR 2 x S and S 1 x B, respectively, and the conformal transformation of M.p onto 
M.p is singular on the bifurcation surface. For this reason the relation between the two 
formulations requires an additional analysis. 

We showed that there exist ultraviolet regularizations which are applicable to both 
covariant Euclidean and canonical free energies. Let us consider for example a situation 
when volume cutoff and Pauli-Villars regularizations are applied simultaneously. Then 
the free energies depend on e and PV parameter //. Because the horizon is excluded there 
is the equality 

F t E [g,{3,^e} = F l c [g,{3^,e} . (5.17) 

As we have shown earlier, the left and right parts of this equality remain finite when e 
cutoff is removed. From (|5.17 ) we obtain 



F?\g,l3,ti) = HmF E [g,p,^e] = limif 0, /i, e] = Ff\g,/3,fi) . (5.18) 

In the limit e — ■> F E becomes the functional defined on M.p, whereas in the same limit 
Ff is defined on A4p. Equation Q5.18Q supports the conclusion Q and |66] that not only 



the divergencies, but also the complete bare free energies in the covariant Euclidean and 
canonical formulations are equivalent when one uses the ultraviolet regularization. 

There are examples where equality (|5.18 ) can be checked by explicit calculations. For 



instance, by rewriting the density of levels (|4.13|) of massive scalar fields in the Rindler 



space-time in the Pauli-Villars regularization one can find from it the canonical free energy 
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and confirm equality ( |5.18| ). An explicit derivation of Eq. (|5.18 ) is also possible in two 
dimensions ||83|| . 

Relation ( |5.18| ) demonstrates that the partition function Z E defined by Euclidean path 
integral ( |2.26|) coincides^ with the canonical partition function Z° , Eq. ( ]2.6|) , including 
the case of spacetimes with horizons. This fact is very important because it justifies the 
statistical-mechanical interpretation of the Gibbons- Hawking path integral method [^8 



59fl when it is applied to black holes. 

It should be noted that ( [5.18 ) was analysed for scalar and spinor fields only The 



relation between the canonical and the covariant Euclidean formulations for high spin 
fields needs an additional consideration. Some discussion of vector fields can be found in 



Refs. 115-117 



5.4 Remarks on rotating and extremal black holes 

So far our discussion was restricted by static black hole geometries. The covariant Eu- 
clidean formulation of statistical mechanics for rotating black holes was studied in Refs. 



118 1 and [94]. It was shown in [ 118 |, an Euclidean manifold which is obtained by the Wick 



rotation of a stationary geometry with the Killing horizon has a conical singularity similar 
to the one which appears in static spaces. This singularity results in the one-loop diver- 
gence which has the form (|4.26|) . The canonical formulation of the statistical mechanics 
of quantum fields on a rotating black hole background was discussed in Refs. [ |119| [-[ |T2"B 1 . 
However, the relation between the canonical and covariant Euclidean formulations in this 
case has not been investigated. 

Some remarks are also in order about extremal black holes. The horizon of an extremal 
black hole has zero surface gravity which indicates that the temperature of such a black 

In the Euclidean theory extremal black holes have the topology 



hole is zero 1241, 125 



of an annulus and there is no conical singularity for an arbitrary period (3 fl 26| [ , [ |T2"7| | . 
The statistical mechanics of quantum fields on the background of the extremal black 
hole has important features which differs it from the non-extremal case. In particular, 



the leading divergence of the density of levels ^ of quantum fields is proportional to 
e L / M where M is the black hole mass and L is the proper distance between the spatial 
boundary and the location of the cutoff near the horizon ||128| |,[ |T2"9|| . Nevertheless, the 
ultraviolet type regularizations still can be applied to eliminate the divergencies of the 
density of states and canonical free energy P9| |. Interestingly, the divergences of ^ and 
F c in these regularizations have the form which is similar to Eqs. fl4.14p -( P:.16| ) and 
( f4.2ip ,( f4\22| ),( |4.2()| ). There are no results comparing the canonical and covariant Euclidean 
formulations of statistical mechanics for extremal black holes. 



6 More precisely, to equate Z and Z E one has to normalize Z E so that to exclude the vacuum energj 
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6 Thermodynamics and statistical mechanics of black 
holes 



6.1 Statistical- mechanical entropy and quantum correction to 
black hole entropy 

Till now we analysed the properties of the statistical-mechanical entropy of quantum fields 
around a black hole. In the general case this entropy is divergent while the Bekenstein- 
Hawking entropy is finite. What is the relation between the statistical mechanical entropy 
of quantum excitations of a black hole and its thermodynamical entropy? 

In quantum field theory the quantum corrections are ultraviolet divergent quantities 
whose divergences are removed by the renormalization of bare coupling constants. As we 
demonstrated the divergencies of the statistical-mechanical entropy S^ v have ultraviolet 
form. It was suggested in []38| that S£ v can be absorbed by the standard renormalition of 



the Newton constant. 

Let us discuss this renormalization approach in more details. The complete informa- 
tion concerning the canonical ensemble of black holes with a given inverse temperature (3 
at the boundary is contained in the partition function Z(/3) given by the Euclidean path 
integral ||59| 

Z(/3) = exp(-I[$]). (6.1) 

Here the integration is taken over all fields including the gravitational one that are given on 
the Euclidean section and are periodic (or antiperiodic) in the imaginary time coordinate r 
with period (3. The quantity $ is understood as the collective variable describing the fields. 
In particular, it contains the gravitational field. Here is the measure of the space 
of fields $ and / is the Euclidean action of the field configuration. The action I includes 
the Euclidean Einstein action. The state of the system is determined by the choice of the 
boundary conditions on the metrics that one integrates over. For the canonical ensemble 
of the gravitational field inside a spherical box of radius r at temperature T one must 
integrate over all the metrics inside tq which are periodically identified in the imaginary 
time direction with period (3 = T~ l . Such a partition function must describe in particular 
the thermal ensemble of black holes. The partition function Z is related to the effective 
action r = — In Z. The free energy F is defined as F = (3~ l T = —f3~ l In Z. 
By using the stationary-phase approximation one gets 

(3F = Y = I[%] + W + .... (6.2) 

Here <3>o is the (generally speaking, complex) solution of classical field equations for action 
/[$] obeying the required periodicity and boundary conditions. Besides the tree- level 
contribution /[$oL expression ( |6.2| ) includes also one-loop corrections W, connected with 
the contributions of the fields perturbations on the background $o ; an d higher order 
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terms in loop expansion, denoted by (...). For free fields, W is a one-loop effective action 
computed in the covariant Euclidean formulation of the theory, see Eqs. ( p. 28 ). 



The one-loop divergences appearing in W can be absorbed by the renormalization of 
the couplings of the initial classical action /. To this aim the latter is chosen in the form 

I(G B ,A B ,c B ) = J d 4 x^L, (6.3) 
+ c B R 2 + CgR^R^ + c B R a i3^ v R al3 ^ v 



(6.4) 

In the presence of one-loop divergences the stationary-phase approximation procedure 
is modified as follows. Denote by W div the UV-divergent part of the one-loop effective 
action W . Then the renormalized quantities are defined as 

I ren = I(G ien , A rcn , clJ = I(G B , A b , 4) + W div , W rcn = W - W dW . (6.5) 

Now the starting point of the semiclassical approximation is in finding the extremal of 
the renormalized action J ren . Since for this background W ian is finite and proportional to 
h, this part of the action describes small quantum corrections. 

The key observation of the renormalization approach is that W div has the same struc- 
ture as (|6.3T) - (|6^f ) and hence W div can be absorbed by simple redefinition of the coupling 
constants in I(Gb, A b ,c b ). In other words, J rcn is identical to the initial classical action 
/ with the only change that the bare coefficients Ab, Gb, and c l B are substituted by their 
renormalized versions A rcn , G ron , and c* cn . W div can be found by using relation (|5.9|) and 
expressions ( |5 . 1 ) - ( 5 . 1 2| ) for the heat kernel coefficients of the corresponding Laplace op- 



erators (on regular backgrounds). The relation between bare and renormalized couplings 
depends on the regularization. For instance, in PV regularization the renormalization of 
the Newton constant for the non-minimally coupled scalar field results is 

1 7^fi-«V • (6-6) 



G Icn G B 2n \6 

where c = In ||| and /i is the PV cutoff. 

The free energy F(/3) = f3~ l Y being expressed in terms of the renormalized constants 
is finite. By calculating F(/3) on the black hole instanton one finds the one-loop free 
energy of a black hole. The "observable" thermodynamic entropy of a black hole S TD [p7 



has the standard form S TD = /3 2 dF(f3)/df3. If we neglect for a while by the logarithmic 
divergencies and put A ron = we obtain 



S IV = S» H (G m ) + 0(h) . (6.7) 

Here S BH (G Icn ) is the Bekenstein-Hawking entropy in the theory of general relativity with 
the Newton constant G ron . The terms 0(h) represent finite quantum corrections propor- 
tional to h. Equations ( |6.6j ) and ( |6.7| ) can be compared with the statistical mechanical 
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entropy S c of quantum fields. For a scalar field the leading divergency of S c can be found 
from (|Og ) 

Sl = ^A ■ (6-8) 

Thus, by using Eq. ( |6.6| ) we find 

S BH (G ien ) = S BH (G B )+S^-Q^ , (6.9) 

where the quantity Q div = (,cfi 2 A/ (2n) appears when a scalar field is non-minimal coupled 
with the curvature. In the next Section we show that Eq. ( |6.9| ) preserves its form when 
not only the leading but all the divergencies are included. 

Equation ( p.7| ) explicitly demonstrates that the "observable" Bekenstein-Hawking en- 
tropy contains the statistical-mechanical entropy of black-hole's quantum excitations as 
its part, but in the general case it does not coincides with it. For non- minimal coupling an 
additional term Q is present. Even in the absence on nonminimal coupling when Q = 0, 
the presence of the bare pure geometrical contribution S BH (G B ) evidently excludes the 
possibility to identify S BH (G Icn ) with S div which has clear statistical mechanical meaning. 
Moreover, in order to have finite value of S BH (G ren ) one needs to assume that pure ge- 
ometrical "entropy" S BH (G B ) is infinite and negative. For this reason the idea to relate 
S BH with quantum excitations does not work, at least in the standard renormalization 
approach. The way out of this problem is to restrict oneself by considering special class of 
the theories where S BH (G B ) = 0. It happens when G" 1 = 0, and hence initially gravity 
is not dynamical. The dynamics of the gravitational field arises as the result of quantum 
effects. The induced gravity is an example of such a theory. 

7 General renormalization and the Noether charge 

Before discussing the models of induced gravity we make comments on the generalization 
of Eq. ( |6.9|) . Let us consider the case when terms quadratic in curvature are preserved 
in the renormalized action. In this case the classical black hole entropy is (see Refs. 
130Hl33j) 



S BH (G Tea , 4J = -±-A - f V^d 2 9 (8nci n R + Airc^Q + Mcljl) . (7.1) 
4G rrn Jt. 



The integral in (|7. 1| ) is taken over the bifurcation surface of the horizon. The first term in 
the r.h.s. of Eq. (|7.1| ) is the Bekenstein-Hawking entropy, and other terms are additions 
because of the high curvature terms in the action. 

For PV and dimensional regularizations in the model with the scalar and spinor fields 
the relation (5J)) takes the formP] 



S BH (G m , 4J = S BH (G B , 4) + Sg - Q div . (7.2) 



7 In Section 4.3 we pointed out that the divergencies S^ iv of the statistical-mechanical entropy in the 
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In PV regularization S div is determined by expression (|4.28|) . The quantity Q div appears 
for nonminimally coupled scalars, 



«- = fsX[ i + "(5- £ ) fl 



(7.3) 



The coefficients a, b depend on PV cutoff \i and are given by Eqs. Q4.23|) , ( }4.24|) . When 
<5div = the general proof of (|7T2|) was given in Refs. 
138H , PSH and |T09| ). 



and 



(see also Refs. [42 



For scalar fields with the nonminimal coupling, quantity ( |7.3|) can be written as 



Qdiv = 



(7.4) 



where it is assumed that the fluctuation of the scalar field ((f) 2 ) is computed in PV regular- 
ization. Relation ( |7.4|) was first found by Solodukhin [f|l]] with the help of the Euclidean 
formulation of the theory with conical singularities. 

The reason why quantity ( |7.4|) appears in formula ( |7.2|) is the following. In the presence 



of scalar fields with nonminimal couplings the Bekenstein-Hawking entropy includes the 



additional term Q = 27r£ J s 2 , where <fr is the classical field, see ||133j| . In quantum theory 



the quantity Q becomes the operator whose average has the divergent part which coincides 
with Q div . In the one-loop approximation Q div determines the quantum correction to S BH 
because of the non-minimal coupling. 



It was shown recently by Wald et al. ||131|| - |l~33| , [ 140|| , ||141|| that the classical black 
hole entropy in diffeomorfism invariant theories of gravity can be interpreted as the 
Noether charge. In Appendix we demonstrate that Q coincides with the Noether charge 
for the nonminimally coupled matter fields propagating on the fixed curved background 
and construct the corresponding Noether current. We also show that 



Q = —(H-E) 

K 



(7.5) 



where k is the surface gravity, H is the canonical Hamiltonian of the fields and E is the 
energy of the fields obtained from the stress-energy tenor, see definitions ( A. 2 ) and ( A. 5 ) 
of Appendix. The latter relation plays an important role in the models of induced gravity. 



volume cutoff method also take the ultraviolet form if one identifies the volume cutoff parameter e with 
an ultraviolet cutoff. The problem is that each field species requires its particular relation between e and 
/x, and there is no universal relation which enables one to remove the divergencies from all fields. For 
instance, f or sca la rs th e relation between the volume cutoff and PV parameters looks as e~ 2 = c^/i 2 , 
see Eqs. (4.11), (4.2S), while for spinors the same relation has to be e~ 2 — c^-fi 2 . This indicates 
that the volume cutoff parameter is not a truly covariant ultraviolet regulator. Thus, regarding the 
renormalization problem, the volume cutoff and "brick-wall" methods run into the difficulty. Some other 

1 137]. For two-dimensional black hole models the 



difficulties of these methods arc discussed in Refs. [134 



relation between the "brick wall" results and the renormalized quantum correction to black hole entropy 
is discussed in Refs. [|S0| and |93||. 
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8 Black hole entropy in induced gravity 



The theory of induced gravity was suggested by Sakharov long ago. The low-energy 
gravitational effective action T[g] in this theory is defined as a quantum average of the 
constituent fields $ propagating in a given external gravitational background g 



exp(-r[s]) = J[D$]exp(-I\g,$]) 



The Sakharov's basic assumption is that the gravity becomes dynamical only as the result 
of quantum effects of the constituent fields. As a result, we have theories of the special 
type, namely theories with I[g] = 0. The gravitons in such a picture are analogous to the 
phonon field describing collective excitations of a crystal lattice in the low-temperature 
limit of the theory. In the general case, each particular constituent field in flS.ip gives 
a divergent contribution to the effective action T[g]. In the one loop approximation the 
divergent terms are local and of the zero order, linear and quadratic in curvature. In 
the induced gravity the constituents obey additional constraints, so that the divergences 
cancel each other. It is also assumed that some of the fields have masses comparable to 
the Planck mass and the constraints are chosen so that the induced cosmological constant 
vanishes. As the result the effective action T[g] is finite and in the low-energy limit it has 
the form of the Einstein-Hilbert action 



V ^\ = -JTr / dVR + 2 dvK\+..., (8.2) 

167rG \Jm JdM J 

where Newton's constant G is determined by the masses of the heavy constituents^. The 
dots in ( |S.2|) indicate possible higher curvature corrections to W[g] which are suppressed 
by the power factors of m~ 2 when the curvature is small. The vacuum Einstein equations 
SW/Sg^ = are equivalent to the requirement that the vacuum expectation values of 
the total stress-energy of the constituents vanishes 

(T^) = 0. (8.3) 



The value of the Einstein-Hilbert action ( |S.2|) calculated on the Gibbons-Hawking instan- 
ton determines the classical free energy of the black hole, and hence gives the Bekenstein- 
Hawking entropy S BH . 

Consider a model of induced gravity |48| that consists of a number of scalar fields 



with masses m s and a number of Dirac fermions with masses m^. Scalar fields can have 
non- minimal couplings £ s . It is convenient to introduce the following two functions 

p(*)=£mf -4$>1*, g(z) = X>*(l-6e.) + 2 5>* . (8.4) 

s d s d 



^ As an another realization of the idea of induced gravity it is worth mentioning the approach by Adler, 



for a review see Refs. [142 ,[L43|. According to Adler, the induced Newton constant may arise as a result 
of the dimensional transmutation because of the dynamical symmetry breaking in the underlying scale 
invariant theory. This picture is different from the mechanism under discussion. 
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Direct calculations show that the induced cosmological constant vanishes when 



p(0)=p(l)=p(2)=j/(2) = . (8.5) 

The induced gravitational coupling constant G is finite if the following constraints are 
satisfied 

g(0)=g(l) = . (8.6) 

Relations (|8.5| ) are satisfied for theories with supersymmetric massive multiplets for which 
p(z) = 0. Equations ( |3.6|) form a linear system defining £ s . 

The presence of the non-minimally coupled constituents is important. In this case it 
is possible to satisfy constraints ( |3.6D on the parameters m s , m<i and £ s which guarantee 
the cancellation of the leading ultraviolet divergencies of the induced gravitational action 
T[g}. The induced Newton's constant in this model is [4S 



\ s d 



•7) 



Let us analyze now the entropy of a Schwarzschild black hole in the induced gravity. If 
constraints ( |3.5| ) and (|8.6|) are satisfied the induced action has logarithmically divergent 
terms only. However these terms are quadratic in curvature and on the Schwarzschild 
background they can be neglected because they do not depend on the geometry. Consider 
the difference S c — Q, where S is the statistical-mechanical entropy of the scalar and 
spinor massive fields (constituents) and Q is the Noether charge which appears because 
of non-minimal couplings of the scalar constituents. S and Q are assumed to be regu- 
larized according to the same prescription, say, by PV method. In the low-energy limit of 
the theory one can calculate these quantities in the leading order in m? . In this approx- 
imation S coincides with the part S^ iv determined in PV-regularization by Eq. (|4.27|) . 
Analogously, the Noether charge is approximated by Q dW , which can be found from Eq. 
( |7.3| ). It is easy to check that for a Schwarzschild black hole the leading divergencies in 
S c — Q are cancelled and from Eqs. ( |4.27| ) and (|7.3| ) one obtains 



S c -Q~S° V -Q** = ^A + C . (8.8) 

Here C is a numerical divergent constant^. In our consideration C can be neglected 
because it does not depend on the black hole geometry. Thus, up to this numerical 
addition, S c — Q agrees precisely with the induced Bekenstein-Hawking entropy S BH = 
jqA. The result Q3.8Q can be also obtained by using the dimensional regularizationQ 

9 The divergence of C is a specific property of the given simplified model which is not free from some 
logarithmic divergencies. In a more complicated model including, for instance, vector massive fields C 
can be made finite. 

10 It should be noted, however, that the check of Eq. (8.8) in the volume cutoff regularization may 



meet a difficulty [144 . Presumably, it happens because this regularization makes the background space 



incomplete, and so using expression (7.4) for Q in the entropy becomes unjustified. 
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As follows from (|8.8|) , the Bekenstein-Hawking entropy is determined by the leading 
order of the statistical-mechanical entropy of the heavy constituent fields. The quantity 
S^iy comes out as the result of integration over a narrow layer located near the horizon 
and having the size of the order of the Comtpton wave length of the constituents, which 
is comparable to the Planck length. Thus, the black hole entropy depends only on the 
local properties of the horizon. This conclusion implies that relation ( |8.8| ) must be also 
valid for a generic static or stationary black holes. 

Another conclusion which can be drawn from ( |8.8|) is that the black hole entropy 
S BH is not identical to the statistical-mechanical entropy S c . Actually, S BH is a finite 
quantity, while S c diverges because each of the constituents gives a positive divergent 
contribution into it. In Eq. ( |8.8|) the divergences of S are eliminated by subtraction of 
the charge Q. As it can be shown, in the model in question Q has positive divergencies. 
Thus, the presence of the fields with non-minimal couplings which is imperative in order 
to make the theory ultraviolet finite also provides the subtraction of the divergencies of 
S c . 

Is there a statistical-mechanical explanation why the Noether charge Q enters relation 
A possible mechanism was suggested in Ref. The induced gravity explains 



the black hole entropy by relating it to statistical-mechanical entropy of the constituents. 
As the result of the backreaction effect, the black hole geometry responds to quantum 
fluctuations of these fields. One consequence of this effect is that fluctuations of the 
energy of fields E near the average value E = cause the fluctuations of the black hole 
mass near its average value M. The black hole entropy S BH determines the number 
density of the states in the interval (M, M + AM) so the quantity exp S BH coincides with 
the degeneracy v(M) of the black hole mass spectrum. Now it is easy to show by using 
differential mass formula H] (see also Eq. (11.2.48) of Ref. [ |145|| ) that in the leading order 



in the Planck constant the change AM of the black hole mass coincides with the change 
AE of the energy of the constituents. Therefore, the problem of finding the degeneracy 
of the black hole mass spectrum can be reduced to the problem of finding the spectrum 
of the energy of the constituent fields, which is more simple. In particular, the black hole 
entropy can be related with the number density v{E) of E as e sBH = v[E = 0). 

Contrary to S BH , the statistical-mechanical entropy S c is related to the spectrum of 
the Hamiltonian H generating translations of the system along the Killing time. H differs 
from E by the Noether charge Q. That is why S BH and S are different and related by 
Eq. (|yD @9|. 



Why the density number of states with the energy E is finite, while the number of 
states of H diverges even in the ultraviolet finite theory? The answer to this question 
is related to the specific property of the quantum system in the presence of the Killing 
horizon. As we pointed out, the spectrum of frequencies of the single-particle excitations 
in this case does not have the mass gap. In other words, there exist modes with negligibly 
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small frequencies, so called soft modes [[15). An arbitrary number of the soft modes can 
be added to such a system without changing its canonical energy H. It is the soft modes 
which result in additional infinite degeneracy of the spectrum of the Hamiltonian H. 
However, contribution of the soft modes to the energy E is not zero and the spectrum of 
E doesn't have the redundant degeneracy 

The reason why the soft modes contribute to E is that E differs from H by the Noether 
charge Q. In fact, Q is determined by the soft modes only. To see this, it is sufficient to 
approximate the black hole geometry (FOJ) near the horizon by the Rindler space (|4.2|) . 
According to Eq. Q7.4Q , Q is determined by the average of the scalar correlator (0 2 ) on the 
bifurcation surface E. In the Rindler approximation the correlator can be easily computed 
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du I d k 



{<p{x)<j){x')) = Tr p<fi(x)(j)(x' 
n ^ 4>u,k{ x )^,k{ x ') + + 1 )0w,fc(a;)Cfe( a;/ ) • ( 8 - 9 ) 



Here n w = (exp(27TC(j//t) — is the number of particles at the Hawking temperature 
with the energy u. The functions 4>u,k{ x ) are the eigen functions of the single-particle 
Hamiltonian H s determined by Eq. Q2.11|) . On the Rindler space 



M x ) = 4^2 (smh(vr W / K )) 1/2 K WK (p(m 2 + h*) 1 ' 2 )e- 4ut e- ik ** j , (8.10) 

where k l and z l are the momentum and the coordinate along the horizon. K iuJ / K are 
the modified Bessel functions. When the proper distance to the horizon p goes to zero 
K iu / K (p{m 2 + A;|) 1//2 ) behave as the delta function 8(u). Thus, when one of the arguments 



of correlator (8.9) is placed on the horizon only the contribution of the soft modes survives. 
By using PV regularization and taking the limit p — > in QS.9|) one gets Q div in form (|7.3|) 
with R = 0. 

The induced gravity proposes the following statistical-mechanical interpretation of 
black hole entropy formula ( |8.8j ). S BH is determined by the density number of physical 
states corresponding to the given black hole mass M. The physical states are related 
to the states of the excitations of the constituents fields in the black hole exterior which 
result in the fluctuation of the black hole mass. The space of physical states is obtained 
by factorizing the space of all thermal excitations over the subspace of the soft modes. 
This removes the additional degeneracy and makes the number density of the physical 
states finite. The factorization over the soft modes is equivalent to subtracting of the 
Noether charge from the entropy S , which reduces the latter to S BH . As was pointed 
out in f4~9 |, there is a similarity between this mechanism and gauge theories, soft modes 



playing the role of the pure gauge degrees of freedom. 
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9 Summary 



This review is devoted to the description of thermal ensembles of quantum fields in a space- 
time of a black hole. This study was stimulated by the attempts to give the statistical- 
mechanical explanation of the Bekenstein-Hawking the entropy. This is the key problem 
of black hole physics. 

The main difficulty of the statistical mechanics of quantum fields near black holes is 
connected with additional thermal (infrared) divergences. In the presence of these diver- 
gences the relation between canonical and covariant Euclidean formulations of the theory 
requires reconsideration. In this review we payed a special attention to the analysis of the 
divergences and methods of their regularizations. It was demonstrated that the canoni- 
cal and covariant Euclidean formulations are equivalent, and in the same regularization 
the divergences in the both formulations are identical. An important property of the 
problem is that the thermal divergences in one regularization take the form of ultraviolet 
divergences in the other regularization. This duality is crucial for the discussion of the 
black hole entropy. Thermal divergences which contribute to the black hole entropy are 
connected with ultraviolet one-loop divergences which renormalize gravitational coupling 
constants. 

We analyzed the problem of black hole entropy and demonstrated it cannot be solved in 
a theory of gravity within the standard scheme of renormalizations. The renormalization 
requires initial bare entropy, which is of pure geometrical origin and (in the absence of 
non-minimal coupling) is infinite negative quantity. However, if the bare classical (tree 
level) gravity is absent the Bekenstein-Hawking entropy S BH can be directly related to 
the statistical-mechanical entropy S of quantum black hole excitations. New important 
feature is that this relation necessarily contains the Noether charge Q for non-minimally 
coupled fields. In one-loop ultraviolet-finite theories without bare gravity there exist the 
relation between S BH and S . 

s bh = s c_q _ ( 91 ) 

These theories belong to the class of so-called induced gravity theories. In these theories 
gravity is induced as the result of collective quantum excitations of heavy constituents of 
the Planckian mass. The same constituent fields which generate low energy gravity are 
responsible for the entropy S BH of a black hole. 

The analysis of concrete models shows that the relation ( |9.1| ) is the result of the 
consistency of the theory. The Noether charge Q determines the difference between the 
energy E of the system and the value of the Hamiltonian H. The entropy S BH describes 
the degeneracy of states of the black hole with respect to its mass. It can be related to the 
degeneracy of the system of constituents with respect to its energy E. The latter differs 
from the degeneracy of the Hamiltonian H. Formula ( |9.1D takes into account this fact. 

In this approach fields (constituents) which contribute into the entropy are to be 
considered as fundamental. One can expect that such fields arise in the fundamental 
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theory of quantum gravity, for instance, in the string theory. This mechanism is not 
known at present, but its possible existence might explain universality of the entropy of 
black holes. 
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A Energy, Hamiltonian and Noether charge 



To explain the meaning of the Noether charge Q discussed in Section 7 let us consider a 
classical real scalar field <ft on the static background with the action 

I[<f>,g] = ~J((f>'^ + m 2 (f) 2 + {;R<f) 2 )V=gd 4 x . (A.l) 

The energy E of the field in the 3D region B is determined by the stress-energy tensor 
T 

E = [ T^da v . (A.2) 

Here da v is the future directed vector of the volume element on B and is the time-like 
Killing vector. The stress-energy tensor is obtained by the variation of action (|A. 1|) 

T„ = ^ . (A .3) 

V^g fig 

In a static spacetime the energy E is conserved on the equations of motion of the field 
(i.e., E does not depend on the choice of B). In addition to the stress-energy tensor 
( |A.3| ), one can define the canonical stress-energy tensor 

dC 

(T c ), u = ^—--g IMl/ C , (A.4) 



>■'■ 

1 1 - 



where C is the Lagrangian of the theory related to the action as / = / y/—gd xC. For 
static spacetimes {T c )^ v yields another conserved quantity which is the Hamiltonian of 
the system 

H = J b T c ^da v . (A.5) 

In the canonical formulation of the theory H plays the role of a generator of the evolution 
of the system along the Killing time. In general the tensors T^ v and (T )^ do not coincide 
and their difference yields the Noether current^ 

j, = -((t%-tJc , (A.6) 

where k is the surface gravity. According to the Noether theorem, this current conserves 
( W M = 0) on the equations of motion. As it follows from ([A. 2 ) and (|A.5| ), the difference 
between the energy E and the Hamiltonian H is the Noether charge Q corresponding to 
the current J M . 

The simplest example when T^ v and {T c )^ u are different is the scalar field with £ ^ 0. 
In this case 

J /t = -^(i?^ 2 + ^(0 2 )^-(0 2 )^)C , (A.7) 

1 1/2 



H-E = U ds k \g 00 \ L/2 (0 2 ), fe - tfw k . (A. 



11 The coefficient 2tt/k in the definition ( |A.6| ) corresponds to the Killing vector ( normalized so that 
C 2 = — 1 at infinity. 
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Here ds k is a three dimensional vector in B normal to the boundary OB and directed 
outward with respect to B. Thus two energies differ by a surface term given on the 
boundary dB of the hypersurface B. Obviously, when one considers a complete Cauchy 
surface the boundary term in (|A.8|) contains only a contribution from the spatial infinity, 
or from the external spatial boundaries if they are present. For a field falling off at 
infinity or obeying suitable conditions at the boundary, E = H. However, the situation 
is qualitatively different if the integration region in E is restricted by the bifurcation 
surface £ of the Killing horizon, where the field can take arbitrary finite values. If 
the contribution from the spatial infinity or external boundary is absent only the Nether 
charge Q on E gives the contribution to the difference H — E 

H-E=^Q , (A.9) 

Z7T 

Q = 2tt£ / V^d 2 6 <f) 2 . (A. 10) 

In an analogous way one can compute the charge Q in other theories. For spin 1/2 fields a 
trivial check shows that the energy and Hamiltonian are identical and Q = 0. Presumably, 
Q is non-trivial in the theories which include fields with higher spins. 
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